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Although the study of ultracold quantum gases trapped by light is a prominent direction of modern 
research, the quantum properties of light were widely neglected in this field. Quantum optics with 
quantum gases closes this gap and addresses phenomena, where the quantum statistical nature of 
both light and ultracold matter play equally important roles. First, light can serve as a quantum 
nondemolition (QND) probe of the quantum dynamics of various ultracold particles from ultracold 
atomic and molecular gases to nanoparticles and nanomechanical systems. Second, due to dynamic 
light-matter entanglement, projective measurement-based preparation of the many-body states is 
possible, where the class of emerging atomic states can be designed via optical geometry. Light 
scattering constitutes such a quantum measurement with controllable measurement back-action. As 
in cavity-based spin squeezing, atom number squeezed and Schrodinger cat states can be prepared. 
Third, trapping atoms inside an optical cavity one creates optical potentials and forces, which are 
not prescribed but quantized and dynamical variables themselves. Ultimately, cavity QED with 
quantum gases requires a self-consistent solution for light and particles, which enriches the picture 
of quantum many-body states of atoms trapped in quantum potentials. This will allow quantum 
simulations of phenomena related to the physics of phonons, polarons, polaritons and other quantum 
quasiparticles. 

PACS numbers: 03.75.Lm, 42.50.-p, 05.30.Jp 



I. INTRODUCTION 

Both quantum optics with nonclassical light and 
physics of ultracold quantum gases are well-established 
and active fields at the forefront of modern quantum 
physicsjl], [2j. However, the integration of the two fields 
is far from being complete. 

Historically, optics treating the light as electromag- 
netic waves is one of the most developed and fruitful fields 
of physics. It has provided a lot of technological break- 
throughs and the highest level of measurement precision. 
Quantum optics extends this to a proper description of 
the quantum fluctuations of light beyond a mean- field de- 
scription. Today it is also a well-established field, both 
theoretically and experimentally [3|. 

The progress in laser cooling techniques in the last 
decades of the 20th century led to the foundation of a 
new field of atom physics: atom optics. It was shown that 
the matter waves of ultracold atoms can be treated sim- 
ilar to light waves in classical optics and can be manipu- 
lated using laser light forces. The quantum properties of 
matter waves beyond the mean-field description became 
apparent after 1995, when the first Bose-Einstein conden- 
sate (BEC) and many other fascination quantum states 
of bosonic and fermionic ultracold atoms were obtained 
[l|, l2j . An exciting demonstration of " quantum atom op- 
tics" was presented in 2002, when the quantum phase 
transition between two states of atoms with nearly the 
same mean density, but radically different quantum fluc- 
tuations was obtained, i.e., the superfluid (SF) to Mott 
insulator (MI) state transition [4|, Q. 



The roles of light and matter in optics and atom optics 
are completely reversed. Various devices known in optics 
as beam-splitters, mirrors, diffraction grating, etc. are 
created using light forces and applied to matter waves. 
However, up to now, the absolute majority of even very 
involved setups and theoretical models in physics of ul- 
tracold quantum gases treat light as an essentially classi- 
cal auxiliary tool to prepare and probe intriguing atomic 
states. In this context, the periodic micropotentials of 
light (optical lattices) play the role of cavities in op- 
tics enabling one to store and manipulate various atomic 
quantum states. 

Quantum optics of quantum gases should close the gap 
between quantum optics and atom optics by addressing 
phenomena, where the quantum natures of both light 
and matter play equally important roles. Experimen- 
tally, such an ultimate quantum level of the light-matter 
interaction became feasible only recently, when the quan- 
tum gas was coupled to the mode of a high-Q cavity 
[g-Q. Even scattering of quantized light from a BEC 
in free (without a cavity) space was not realized so far. 
However, implementing cavity quantum electrodynamics 
(QED) setups with quantum gases as recently achieved 
provides for the most interesting and controllable imple- 
mentation to study the atom- and light-generated quan- 
tum effects. 

On the one hand, the quantum properties of atoms will 
manifest themselves in the scattered light, which will lead 
to novel nondestructive methods of probing and manipu- 
lating atomic states by light measurement. On the other 
hand, the quantization of light (i.e. trapping potentials) 




FIG. 1: Setup. The atoms periodically trapped in a lattice 
are illuminated by the transverse probe ao (the trapping laser 
beams are not shown). The scattered light mode a\ is col- 
lected by a cavity. Another probe r/i through a cavity mirror 
can be present. The photons leaking the cavity are measured 
by the detector. 



will modify atomic many-body dynamics well-known only 
for classical potentials and give rise to novel quantum 
phases. 

This paper is organized as follows. In Sec. II, we for- 
mulate a general theoretical model, which treats the light 
and atomic motion in a fully quantum way. As a partic- 
ular example, a generalized Bose-Hubbard model will be 
derived for quantized optical lattice potentials. Section 
III considers quantum non-demolition (QND) setups for 
probing the quantum states of ultracold particles by light 
scattering. In this section, we establish the relations be- 
tween various light- and matter-related observables (such 
as expectation values of different quantities, their distri- 
bution functions and correlation functions). In Sec. IV, 
we go beyond the standard goal of the determination of 
correlation and distribution functions: we consider the 
truly quantum nature of the measurement process itself, 
which is evident in a single run of an optical measurement 
without averaging over many realizations. The quantum 
measurement back-action is used as an active tool to pre- 
pare various many-body states of ultracold particles. In 
Sec. V, we proceed deeper into the quantum regime 
and analyze the situation, where not only the quanti- 
zation of probe light and atomic motion are important, 
but even the quantum nature of the trapping potential is 
crucial. This brings us to the formulation of the concept 
for "quantum optical lattices" produced inside a high- 
Q cavity. Various applications of the cavity QED with 
quantum gases are mentioned in Sec. VI, and conclusions 
are presented in Sec. VII. 



II. THEORETICAL MODEL 

Here, we formulate the theoretical model to describe 
the interaction of an ultracold quantum gas with quan- 
tized light |9l-[l2| • The spatial geometry of light is taken 
into account in a very general form: the model can de- 
scribe the interaction inside a cavity or in free space. The 



trapping of atoms by fully quantum potential is captured 
by this model as well. Describing the quantum properties 
of atoms, we focus on a particular case: spinless bosons. 
Generalizations to fermions and spin particles is straight- 
forward and of interest as well. 



We consider an ensemble of N two- level atoms in an op- 
tical lattice with M sites. In quantum optics of quantum 
gases, a restriction to optical lattices is very convenient 
as it allows to precisely describe the many-body atomic 
state for a broad range of parameters. Moreover, exper- 
imentally different setups can be described by the same 
general model. The typical examples include a multi-site 
lattice with the low filling factor (e.g. one or two atoms 
per lattice site as in a typical superfluid to Mott insulator 
transition setup) and a trapping potential with the high 
filling factors but small number of sites (e.g. a BEC in a 
double- well potential setup). 



In general, the atoms are trapped in a lattice potential 
created by strong lasers as it is usual in the standard cold 
atom problems. In the presence of such a classical poten- 
tial, atoms are illuminated and scatter light at different 
directions. As shown in Fig. [TJ the atoms in a lattice 
(green) are illuminated by a probe beam (red), and the 
measurements are carried out in the direction of one of 
a scattered light modes (blue). In fact, the probe and 
scattered light modes can be in free space without the 
presence of any cavity. In practice, light modes can be 
selected by traveling- or standing-wave cavities, or even 
correspond to different modes of the same cavity. One 
important reason for including a cavity is to enhance 
light scattering into some particular direction. Another 
reason, is that the cavity mode can form a fully quan- 
tum trapping potential for the atoms (even if no classical 
trapping potential is present). For defmiteness, we will 
consider the case, where the light mode functions are de- 
termined by cavities, whose axes directions can be varied 
with respect to the lattice axis (the simplest case of two 
waves, probe beam and cavity mode, at angles 9 and 6\ 
is shown in Fig. [1]). Instead of varying the angles, also 
the light wavelengths can be varied with respect to the 
wavelength of a trapping beam. We also assume, that 
not all M lattice sites are necessarily illuminated, but 
only a subset of K < M sites. The selection of K out 
of the total M sites enriches the picture. In the simplest 
case, a continuous part of a lattice with K sites can be 
illuminated. However, the nontrivial selection of the illu- 
minated sites is also possible: e.g., each second site can be 
easily illuminated by choosing the light wavelength twice 
as the wavelength of the trapping beam. Moreover, using 
several probe beams one can make the optical geometry 
even more interesting. 



A quite general model many-body Hamiltonian in the 



second quantized form is given by |Tol - [T2| 

H = Hf + H a , 



2Tta s h 



H a = /"dW(r)iI al *(r) 



d 3 r&(r)&(r)V(r)V(r). 



(la) 
(lb) 



(lc) 



In the field part of the Hamiltonian Hf, ai are the an- 
nihilation operators of light modes with the frequencies 
loi , wave vectors k; , and mode functions ui (r) , which can 
be pumped by coherent fields with amplitudes r\i. In 
the atom part, H a , ^>(r) is the atomic matter-field op- 
erator, a s is the s-wave scattering length characterizing 
the direct interatomic interaction, and H a i is the atomic 
part of the single-particle Hamiltonian H\ , which in the 
rotating-wave and dipole approximation has a form 



H 1 =H f +H al ,(2&) 



huj a 



H a \ = h -7r<7z ~ ihy^[a + giaiui(r) - H. c.](2b) 



Here, p and r are the momentum and position opera- 
tors of an atom of mass m a and resonance frequency u a , 
<j + , cr~ , and a z are the raising, lowering, and population 
difference operators, gi are the atom-light coupling con- 
stants for each mode. The inclusion of the interaction be- 
tween the atom and quantum light in the single-particle 
Hamiltonian is the key step, which is different from the 
standard problems of ultracold atoms in classical poten- 
tials. 

We will consider essentially nonresonant interaction 
where the light-atom detunings A; a = u>i — uj a are much 
larger than the spontaneous emission rate and Rabi fre- 
quencies giai. Thus, in the Heisenberg equations ob- 
tained from the single-atom Hamiltonian Hi ©, the 
atomic population difference a z can be set to —1 (ap- 
proximation of linear dipoles, i.e., the dipoles responding 
linearly to the light amplitude with the negligible popu- 
lation of the excited state). Moreover, the polarization 
<j~ can be adiabatically eliminated and expressed via the 
fields ai. The elimination of er~ and setting a z = — 1 
is equivalent to the adiabatic elimination of the upper 
atomic level. The light-atom detunings can be all then 
replaced by a single value A a = oj p — u> a , where uj p is, 
for example, the frequency of the external probe. An 
effective single-particle Hamiltonian that gives the cor- 
responding Heisenberg equation for a; can be written as 
-Hicff = H f + H a i with 



H a i 



2m a 



V d (r) 



A a 



^u* (r)u m (r)gig m aja m . (3) 



Here, we have also added a classical trapping potential of 
the lattice, V^i(r), which corresponds to a strong classical 



standing wave. This potential can be, of course, derived 
from one of the modes ai = a c \ [in this case Ki(r) — 
hq^\a r ,]U r ,](r)\ 2 / A r ]„], and it can scatter light into other 
modes J13J. Nevertheless, at this point we will consider 
V c \(r) as an independent potential, which does not affect 
light scattering of other modes that will be significantly 
detuned from a c \ [i.e. the interference terms between a c \ 
and other modes are not considered in the last term of 
Eq. (]3j)]. A later inclusion of the light scattered by the 
trapping wave will not constitute a difficulty, due to the 
linearity of dipoles assumed in this model. 

If one considers scattering of weak modes from the 
atoms in a deep lattice, the fields ai are much weaker 
than the field forming the potential V^(r). To derive the 
generalized Bose-Hubbard Hamiltonian near zero tem- 
perature, we expand the field operator ^(r) in Eq. (p}, 
using localized Wannier functions corresponding to K;i(r) 
and keeping only the lowest vibrational state at each site: 
\p(r) = 5Z i=1 biiu(r — Ti), where bi is the annihilation op- 
erator of an atom at the site i with the coordinate r^. 
Substituting this expansion in Eq. ([TJ with H a i ([3]), we 
get 

M 

H = Hf+Y J JtA h i 



I. . m 



V '-J = l 



u 



f;6t6i(6j6i-l), 



(4) 



where the coefficients J^ correspond to the quantum mo- 
tion of atoms in the classical potential and are typical for 
the Bose-Hubbard Hamiltonian [4(: 



./: 



'•j 



drw(r 



h 2 V 2 

2m 



Kl(r) w(r 



(5) 



However, in contrast to the usual Bose-Hubbard model, 
one has new terms depending on the coefficients J^J 1 , 
which describe an additional contribution arising from 
the presence of light modes: 



•/. 



lm 



drw( 



)u* (r)u m (r)w(r 



(G) 



In the last term of Eq. (j4]), only the on-site interaction 
was taken into account and U = 47ra s 7i /m a J dr\w(r)\ 4 . 
Note, that if the contribution of the quantized light 
is not much weaker than the contribution of the classi- 
cal potential, or if the classical potential is not present 
at all, the Wannier functions should be determined in a 
self-consistent way taking into account the mean depth of 
the quantum potential, generated by the quantum light 
modes. This indeed significantly complicates the theo- 
retical treatment [lSq . Although the form of the above 



equations still holds, the coefficients ([5]) and ([6]) will de- 
pend on the quantum state of light as well. 

As a usual approximation, we restrict atom tunneling 
to the nearest neighbor sites. Thus, coefficients ([5]) do 
not depend on the site indices (Jf\ = J§ x and Jf i±1 = 
J cl ), while coefficients © are still index-dependent. The 
Hamiltonian ((4]) then reads 



H = H f + J C W + J C[ B 

/ K \ 



sr 



Lm 



t I V^ Tin 

I Urn I 2_^ J i,i 



vi=l 



K 



j-^9i9ma\a m ^ J iJ b t b J 



M 

where < i,j > denotes the sum over neighboring pairs, 
hi = bjbi is the atom number operator at the i-th 
site, and B = J2i=i &J&*+i + H.c. While the total atom 
number determined by N — 5Z i=1 ^ is conserved, the 
atom number at the illuminated sites, determined by 
Nk = X)j=i ^i: i s n °t necessarily a conserved quantity. 

The Heisenberg equations for ai and bi can be obtained 
from the Hamiltonian ([7]) as 



ai = —i I wi + 



ML 

A a 



A" 



K 



E^+f E^ N 



i=l 



A a 



<M> 
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_i X" 7 . SmOm | 2^ J i,T^ 



m^l 



■ 9i 
''A a 



K 



^2g m a m ^ J i7 b i b i +»7; - K/a;,(8a) 



m^l 



K <i,j> 



k = -- Jq + -7-^2gig m <4amJi™ + Uhi I 6,; 



I :>n 



* / tcI "" V^ t Tim I 
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where we phenomenologically included the decay rate 
Ki of the mode a;. We do not add the corresponding 
Langevin noise term, since we will be interested in the 
normal-ordered quantities only. The decay of the atoms 
can be included in a similar way as well, but it is usually 
much smaller than the cavity relaxation. 

In Eq. (|5aj) for the electromagnetic fields a;, the two 
last terms in the parentheses in the first line correspond 



to the phase shift of the light mode due to nonresonant 
dispersion (the second term) and due to tunneling to 
neighboring sites (the third one). The second line in 
Eq. (|5aj) describes scattering of all modes into ai , while 
the third term takes into account corrections to such scat- 
tering associated with tunneling due to the presence of 
additional light fields. In Eq. f|8b[) for the matter field 
operators bi , the first line gives the phase of the matter- 
field at the site i, the second and third terms describe 
the coupling to neighboring sites. 

It is important to underline that except for the di- 
rect coupling between neighboring sites, as usual for the 
standard Bose-Hubbard model, Eqs. ((8]) also take into 
account long-range correlations between sites, which do 
not decrease with the distance and are provided by the 
common light modes a; that are determined by the whole 
set of matter-field operators bi. Such a cavity- mediated 
long-range interaction and nonlocal correlations between 
the operators bi, which are introduced by the general 
Eqs. (O), can give rise to new many-body effects beyond 
the standard Bose-Hubbard model 1121. 



III. QND PROBING OF THE MANY-BODY 
STATES OBSERVING LIGHT 



As the results of the previous section show, the light- 
matter interaction leads to the joint evolution of the light 
and atomic variables. The light and matter get corre- 
lated. As this is a fully quantum problem, in general, 
the light and matter get entangled. Thus, observing the 
light, one can obtain the information about the quantum 
states of the atoms. The goal of this section is to es- 
tablish relations between the characteristics of scattered 
light and those of the many-body atomic system. 

Considering measurements in quantum mechanics, a 
first goal is to obtain the expectation values of desired 
quantities or their distribution functions. Thus, such 
type of probing is intrinsically associated with multi- 
ple measurements of the same quantity and averaging 
over the measurement outcomes. The type of the quan- 
tum measurement considered in this section corresponds 
to the quantum nondemolition (QND) [ljj observation 
of different atomic variables by detecting light. In con- 
trast to the completely destructive schemes usually used 
(e.g. time of flight measurements, absorption images), 
the QND measurement by light scattering only weakly 
perturbs the atomic sample and many consecutive mea- 
surements can be carried out with the same atoms with- 
out preparing a new sample for each measurement. 

However, an important statement of quantum mechan- 
ics is, that any measurement, even a QND one, affects 
the quantum state of the system. Therefore, to measure 
the expectation values or statistics of some variable, one 
should prepare the system in the same quantum state 
before each measurement (or wait until the initial state 
will restor due to the system evolution). In the next sec- 
tion, we focus on the essentially quantum properties of 



the measurement process itself going beyond the simple 
goal of measuring only the expectation values. We will 
analyze, how the many-body atomic state changes during 
the measurement process, i.e., we will consider the mea- 
surement at a single quantum trajectory (single run of 
a measurement) without a requirement of the statistical 
averaging over the many runs of an optical experiment. 
Let us first consider the model of spinless bosons intro- 
duced before to demonstrate, how the relation between 
the light and atom observables can be established. Then, 
we will provide the generalizations of that systems and 
measurements suggested recently. None of those exam- 
ples goes beyond the typical measurements of expectation 
values and statistics utilizing the averaging procedure. 
The underlying quantum properties of the measurement 
process (measurement back-action at a single quantum 
trajectory) will be addressed in the next section. 



A. Model for probing bosons in optical lattices 

To focus on the question how to establish the relation 
between the light and atom observables, we simplify the 
model. We consider only two light modes: the probe 
a and the scattered light a\. We assume that the tun- 
neling of atoms between the neighboring sites is much 
slower than light scattering, and tunneling can be ne- 
glected during light scattering. Physically this means 
that the quantum properties of the atomic state are de- 
termined by tunneling and interaction, but are frozen 
during the measurement time. In practice, the tunneling 
and scattering times can be different in orders of mag- 
nitude. For a deep lattice the coefficients J\™ reduce to 
J-™ = u ; *(r i )u m (r i ) neglecting atom spreading. The in- 
fluence of the atom spreading within each site on the light 
signal can be studied even by classical scattering [l5[ . 

After those simplifications, the Hamiltonian Eqs. (|7|) 
takes the form: 

H = Ti(lui + UiiDii)a[ai + hU w (D* w a* ai + D w a a\) 

-ihijilai - 77iaJ),(9) 

where a\ is the cavity-mode annihilation operator. The 
external probe is assumed to be in a coherent state, thus 
its amplitude is given by a c- number a,Q. £/j m = gig m /A a 
(Z, 77i = 0,1), 771 is the amplitude of the additional prob- 
ing through a mirror at the frequency uj p (the probe- 
cavity detuning is A p = u p — Ui). The operators 
Dim — Y]j—i u i{ r j) u m(Yj)nj sum contributions from all 
illuminated sites with the atom-number operators hj at 
the position Tj. 

The first term in the Hamiltonian describes the atom- 
induced shift of the cavity resonance. The second one re- 
flects scattering (diffraction) of the probe ao into a cavity 
mode a\. The key feature of a quantum gas is that the 
frequency shift and probe-cavity coupling coefficient are 
operators, which leads to different light scatterin g am pli- 
tudes for various atomic quantum states [l(J, LLll, LLil LL7| ■ 



The Hamiltonian ([9]) describes QND measurements of 
the variables associated with the operators L>i m by de- 
tecting the photon number a\ax or light amplitude a\ 
related quantities (e.g. quadratures). In order for a mea- 
surement to be a QND one, several conditions should be 
fulfilled for the "signal observable" of interest As {Di m 
in this case), "probe observable" Ap, which is actually 
detected, (here, a\a\ or a\) and the coupling between 
them through the interaction Hamiltonian [14J . Accord- 
ing, for example, to Ref. [lj], the interaction Hamilto- 
nian should be a function of As ; the interaction between 
the signal and probe should affect the dynamics of Ap 
([Ap,H] 7^ 0), whereas the signal observable should not 
be affected by the coupling to the probe ([As,H] = 0). 
In addition, measuring As, its conjugate variable is af- 
fected in an uncontrollable way, therefore the evolution 
should not depend on that uncontrolled variable at all. 
In this case, the signal observables L>i m depend only on 
the atom numbers, and their conjugate variables related 
to the atomic phase are not present in the Hamiltonian. 
One can see that all those conditions of the QND mea- 
surement are fulfilled for the Hamiltonian ©. 

Note, that one has a QND access to various many- 
body variables, as Di m strongly depends on the lattice 
and light geometry via uo.i(r). This is an advantage 
of the lattice comparing to single- or double-well setups 

For example, one can consider a ID lattice of the pe- 
riod d with atoms trapped at Xj — jd (j = 1,2, ..,M). 
In this case, the geometric mode functions can be ex- 
pressed as follows: uo,i(r,-) = exp(ijko t i x d + 4>o,ij) for 
traveling waves, and uo,i(rj) = cos(jko,xxd + (f>o,ij) for 
standing waves, where ko t xx = |ko,i| sfn#o,i) $0,1 are the 
angles between mode wave vectors ko,i and a vector nor- 
mal to the lattice axis. In the plane- wave approximation, 
additional phases 0o,ij are j-independent. 

For some geometries, Dxx simply reduces to the op- 
erator Nk = X^7=i^j °f the atom number at K sites 
[lOl [ill , [la ] (if ai is a traveling wave at an arbitrary 
angle to the lattice, or the standing wave with atoms 
trapped at the antinodes). If the probe and cavity modes 
are coupled at a diffraction maximum (Bragg angle), 
i.e., all atoms scatter light in phase, u*(tj)uo(vj) = 1, 
the probe-cavity coupling is maximized, Dxo — Nk- 
If they are coupled at a diffraction minimum, i.e., the 
neighboring atoms scatter with opposite phases and n, 

Dxo = Z),- =1 (-l) J ' +1 ftj = ^odd-^Veven is the operator of 
number difference between odd and even sites. Thus, the 
atom number as well as number difference can be non- 
destructively measured. Note, that those are just two of 
many examples of how a QND- variable can be chosen by 
the geometry in a many-body system. 

From the Hamiltonian, the Heisenberg equation for the 
scattered light can be obtained as follows: 



«i 



■i (u 



u\ + UiiDn ax - iUxoDxoa - nax + r?i,(10) 



where n is the cavity decay rate. In classical physics, 
such an equation is directly analogous to the Maxwell's 
equation for the light amplitude of the cavity mode, and 
the classical meaning of the Di m operators (i.e. the fre- 
quency dispersion shift and coupling coefficient between 
two modes) is obvious. Here, in the fully quantum prob- 
lem, both light- and matter-related quantities are treated 
as operators indeed. 

The stationary solution for the operator of the light 
amplitude oscillating at the probe frequency takes the 
form 



the mean atom numbers only. However, the number of 
photons scattered into a cavity depends quadratically on 
the atom operators: 



a\ = 



Vi - UJioOoDio 
i(Unt>u - A p ) + 



(11) 



which gives us a direct relation between the light operator 
and various atom number-related operators. It is clear 
that if a light-related observable is a linear function of the 
atom number operators hj , then the measurement of that 
observable will depend only on the mean atom numbers 
(i.e., their expectation values {hj)). Such a measurement 
will carry information only about the mean atom density, 
which can be similar for various quantum states, and thus 
is not of interest for the scope of this paper. Therefore, 
the question is to find light observables, which depend 
nonlincarly on the atom number operators. In this case, 
the measurement will reveal the higher moments of the 
atom number operator, which carry information about 
the quantum state of ultracold atoms. 

This suggests us to consider the following optical con- 
figurations, where the light observables are sensitive to 
various atomic states. 

(I) Transverse probing [1(J, [ll], [ItJ . Here we neglect the 
dispersive frequency shift. In this case, the light ampli- 
tude (ai) is a linear function of the atom numbers and is 
not of great interest. However, the light intensity (given 
by the mean photon number n$ = {a\a\)) already de- 
pends on the atom density-density correlations {hihj), 
which differ for various atom states. Moreover, the pho- 
ton number variance carries the information about the 
four-point correlation function (hihjhkhi) , which is even 
a more exciting result. 

(II) The probing through a mirror, where the disper- 
sive frequency shift plays a key role [la . Il7j . In this 
case, even the light amplitude nonlinearly depends on the 
atom numbers (see the denominator in Eq. (fTTj) ). Here, 
the measurement of light can directly reveal the full atom 
number distribution function. The measurement of fre- 
quency shifts in a cavity configuration is similar to the 
measurements of the phase shifts in the free-space geom- 
etry. 

Let us consider the configuration (I), where only the 
transverse probe ao is present (no probing through a mir- 
ror 7/1 = 0) and the dispersive frequency shift is small (the 
term UuDu can be neglected) in Eq. (jlll) . The light am- 
plitude operator is linear in the atom number operators, 
ai = CDio = J2j=i Ajfij, where C = iUioao/(iA p — n) 
and Aj — u*(rj)uo( r j)- Its expectation value measures 



n$ 



(a\ai) 



K 

\c\ 2 {b* w b w ) = \c\ 2 J2 AiA.ih^xn) 



Importantly, already a simple quantity as the mean light 
intensity (mean photon number n.$) depends on the atom 
density-density correlations, which are second moments 
of the atom operators and arc different for various atomic 
quantum states. 

Similarly to the light amplitude, the quadratures of 
the light fields are linear in atom numbers. However, the 
quadrature variances depend on the density-density cor- 
relations and thus are sensitive to the atomic state. As 
the mean photon number is very sensitive to the second 
moments of the atom numbers, the photon number vari- 
ance (given by the second moments of n$) depends on the 
fourth moments of the atomic operators {hihjhkhi) and 
provides even more information about the atomic state 
through the four-point correlations. Thus, the light in- 
tensity (number of scattered photons), quadrature vari- 
ance, and photon number variance all are sensitive to 
various atomic states with different atom number fluctua- 
tions. Here we focus only on the intensity measurements, 
while the quadratures and photon number variance were 
considered in details in |ll| . 

For the simplest case of two traveling light waves, the 
situation is almost identical to the simple diffraction in 
free space: the probe ao is scattered from the periodic 
grating of atoms and the outgoing light a\ is measured 
at some angle to the lattice. Consequently, the angular 
distribution of the scattered light amplitude (ai) shows 
the classical diffraction pattern. In the diffraction max- 
ima (Bragg angles) , all atoms scatter light in phase with 
each other (all Aj — 1) and the light amplitude is pro- 
portional to the mean atom number in the illuminated 
region Nk- In the diffraction minima, the mean light 
amplitude is zero due to the destructive interference (the 
neighboring atoms scatter light with the opposite phases 

A, = (-iy +1 ). 

However, the photon number exhibits a more interest- 
ing angular distribution. For a quantum gas, the correla- 
tion function {hihj) generally cannot be factorized. As a 
consequence, the photon number is not given by the light 
amplitude squared, but is an essentially different quan- 
tity. For example, in the diffraction maximum the pho- 
ton number is proportional to {N^) ^ Nf^. Furthermore, 
the light intensity in the diffraction minimum is not zero 
at all, but reflects the quantity ((./Vodd — N cvcn ) 2 ) ^ 0. 
The latter quantity corresponds to the fluctuations of the 
atom number difference between the odd and even sites. 
Although the atom number difference is zero in average 
(the light amplitude (ai) is zero as well), its fluctuations 
are very different for various atomic states, which is re- 
flected in the mean photon number n$. 

To characterize the difference between the quantum 
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FIG. 2: Angular distribution of scattered light. Two traveling waves are used for probing and measurement, (a) Intensity of 
classical scattering shows usual diffraction peaks, (b) Quantum addition to classical scattering R(8i). While the Mott insulator 
(MI) state shows no (zero) addition, the quantum addition for scattering from the superfluid (SF) state shows the anisotropic 
signal proportional to the atom number N. N = M = K = 30, lattice period is d — A/2, the probe angle is 9q — 0. 



and classical contributions to the light scattering, we in- 
troduce a quantity proportional to the difference between 
the photon number and the amplitude squared: 

R(e ,e 1 ) = (b*b)-\(b)\ 2 = 

K 

= (6n a Sn b )\A\ 2 + ((Sn 2 ) - (6njn b )) J^ \ A if, (13) 



which is the angle-dependent function of various atom 
number fluctuations. Here, A(6q,6i) — J2i=i -^(^cb #i) 
is the angular distribution of classical diffraction. 

The atom number fluctuations are the following: Sfii = 
fii — n giving the correlations between different sites 
(Sh a 5fib) — (fi a fib) —n 2 , and the on-site fluctuations (Sh 2 ) 



equal to the variance (An,) 



Although the 



approximation that pair fluctuations are equal for all sites 
is not general, it hold for several interesting examples 
such as Mott insulator state (MI), superfluid state (SF) 
and the coherent state approximation to the SF state. 
Importantly, all those states have the same mean atom 
numbers n at each site, but very different atom num- 
ber fluctuations. We will show that such differences are 
readily captured by the light scattering. 

In the MI state, there are no atom number fluctuations, 
and the light intensity shows precisely classical diffraction 
pattern, thus the quantum - classical difference is zero, 
^?mi = 0. In the approximate coherent state, the pair cor- 
relations between the sites are neglected, (8n a Sn b ) = 0, 
but the on-site fluctuations are Poissonian, (5h 2 ) = n. 
Hence, in Eq. (|13l) the first term with classical angular 

distribution is zero and Rcoh — n X)i=i l^il 2 - Note that 
the angular distribution here is completely different from 
the classical diffraction pattern, as the sum does not go 
over the light amplitudes, but rather over light ampli- 



tudes squared (i.e. the intensities of laser beams at the 
positions of the atoms). In the SF state, both on-site 
and pair fluctuations are non-zero and both two terms 
in Eq. (|13p with two different angular distributions con- 
tribute to i?SF- 

Figures 2 shows the angular distributions of two quan- 
tities: the light amplitude and quantum - classical dif- 
ference ("quantum addition" for the classical signal) R 
for the SF state. In this figure, a case of two traveling 
waves is presented. In classical scattering, two diffrac- 
tion maxima are present, while outside them, there is a 
region of the destructive interference. However, the angu- 
lar dependance of R is different: it is nonzero at diffrac- 
tion minima and zero at the directions of the diffraction 
maxima. Thus, the photon number is nonzero every- 
where. In the diffraction maxima, it is proportional to 
N 2 , while outside them it is isotropic and proportional 
to the atom number variance, which for the SF state is 
N. The scattering from the MI state gives zero photon 
number outside the diffraction maxima. 

One should underline an important property of light 
scattering from periodic lattices, which is very convenient 
for experiments. Although the photon numbers in max- 
ima and minima are different in the order in atom num- 
ber, there is no problem to resolve both of them, because 
they are completely separated by the angles. The quan- 
tum contribution proportional to N appears on top of a 
zero classical contribution. In the case of scattering from 
a homogeneous BEC, this contribution would appear on 
the top of the strong classical forward scattering, propor- 
tional to A^ 2 , and could be hardly resolved. 

A coherent product state of atoms in different wells is 
an approximate description of the SF state, where pair 
correlations are neglected and, as a consequence, the to- 
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tal atom number fluctuates. Although the latter fact 
is a definite defect of the theory, such a coherent-state 
approximation is usually used in the many-body theo- 
ries, for example, as a step to derive the mean-field ap- 
proach. It is important, that light scattering is sensitive 
to such an approximation and can distinguish between 
the SF and coherent product state. In Fig. 2, the func- 
tion R would show absolutely isotropic angular distribu- 
tion without going to zero at the Bragg angles. Indeed, 
Eq. (fT3"|) demonstrates that R is sensitive to the pair cor- 
relations and the total atom number conservation. 

Figure 3 shows the same quantities as in Fig. 2, but 
for the configuration, where two standing waves are used 
instead of two traveling ones. The classical scattering 
consists of four diffraction maxima, which can be simply 
explained by the fact that each of two standing waves 
effectively consists of two traveling waves. However, the 
striking feature of the quantity R is that its angular dis- 
tribution is much richer than the classical scattering. The 
scattering from a quantum gas displays features at the 
angles, where the classical scattering does not exist at 
all. More precisely, in classical scattering only the zero- 
order diffraction maxima are possible for the parameter 
set considered. However, scattering from a quantum gas 
displays the features at angles, which would correspond 
to the first-order diffraction maxima (i.e. they appear at 
the angles corresponding to the doubled lattice period) . 

Such a behavior can be explained by the last term in 
Eq. ([TBf , which sums the light intensities at the positions 
of the atoms. For the traveling waves, the mode func- 
tions are exponents and all \Ai\ 2 = 1 giving the isotropic 
contribution. However, for the standing waves, the mode 
functions are cosines, and their squares display the period 
doubling and continue to depend on the position of the 
atoms. In Eq. (fTB"]) . this geometrical term is multiplied 
by the on-site and pair fluctuations. Physically, one can 
say that scattering from or into a standing wave reflects 
not only the periodic density, but also the periodicity of 
the density fluctuations. 

More generally, one can show that the multiple features 
in the fluctuation-dependent scattering appears not only 
for the standing waves, but for any scattering which in- 
volves the amplitude modulation inside the sample. As 
shown in [ll| , similar angular features can be obtain if the 
quadratures are observed (even for two traveling waves) , 
because different quadratures are coupled spatially dif- 
ferently to the atoms in a lattice. The angular distri- 
bution of photon statistics displays similar features as 
well. Another way to observe those features, and thus to 
probe the periodic structure of the atomic noise, is to use 
several probe beams, which can be traveling waves. In 
this case, several beams will interfere inside the sample, 
which results in the amplitude modulation within the lat- 
tice. The latter situation can, in principle, be realized in 
a recent experiment [3lJ on light scattering from the ul- 
tracold atoms in optical lattices, where five probe beams 
were used to simultaneously cool and probe the atoms. 
Together with the work on light scattering from atoms 



in a 3D lattice [l5j, these are the first experiments on 
light scattering from ultracold atoms in optical lattices 
in there truly quantum regime, where the main object of 
interest and measurement is the light, rather than atoms 
in contrast to most of the works involving time-of-flight- 
type measurements. 

Let us now turn to the configuration (II) mentioned 
above: probing through a mirror without any transverse 
probe. In Eq. (jlll) . we set ciq = 0, and the photon 
number operator a\a\ will be therefore represented by 
a lorcntzian function with a frequency shift given by the 
atom number-dependent operator U\\D\\. In the simple 
cases (the traveling-wave cavity or standing wave cav- 
ity with atoms trapped at their antinodes), the operator 
D\i reduces to the operator of the atom number inside 
the cavity Nk- This type of measurement corresponds 
to the detection of cavity transmission spectrum. For a 
given detuning between the probe and cavity frequencies 
A p , the expectation value of the photon number is 



n$(A p ) 



OJOl) 



\m\ 



(A p - U 1X N K ) 2 



(14) 



where the expectation value is calculated for the many- 
body atomic state. This expression shows that the pho- 
ton number depends nonlinearly on the atom number 
operator and, thus, is sensitive to the quantum state of 
atoms. 

In the atomic state, where the atom number in K sites 
does not fluctuate and takes only a single value Nk, the 
above expression is simply a classical Lorentz contour 
describing the cavity transmission. This is the case for 
the MI state. However, in the states, where the atom 
number fluctuates, the cavity transmission is much more 
complicated. As was shown in 16] for the good cav- 
ity limit, the transmission spectrum consists of multiple 
lorentzians forming a comb-like structure (see Fig. 4 for 
the comparison of MI and SF states). The peaks corre- 
spond to all possible atom numbers in a cavity, thus, the 
peaks appear at all possible dispersion frequency shifts. 
Moreover, as was proved in [la ], such a transmission 
spectrum directly maps out the full atom number dis- 
tribution function. Thus, by measuring the transmission 
spectrum of a cavity containing ultracold atoms in some 
particular quantum state, one can directly obtain the full 
atom number distribution function corresponding to that 
quantum state. 

Figure 4 shows the cavity transmission spectra (and 
hence directly the form of the atom number distribution) 
for MI and SF states. The atom number fluctuations are 
completely suppressed in MI state and one sees only a 
single classical Lorentz contour. In SF state, the atom 
number distribution is close to the Poissonian one (and 
hence close to the Gaussian one for some conditions) . In 
Fig. 4(c), the transmission spectra are shown for differ- 
ent number of sites K illuminated. One observes how 
the transmission spectrum and the atom number distri- 




FIG. 3: Angular distribution of scattered light. Same as in Fig. 2, but two standing waves are used for probing and measurement, 
(a) Intensity of classical scattering shows usual diffraction peaks, (b) Quantum addition to classical scattering R(0i). While 
the Mott insulator (MI) state shows no (zero) addition, the quantum addition for scattering from the superfluid (SF) state 
shows the anisotropic signal proportional to the atom number N. In contrast to the traveling wave configuration (Fig. 2), the 
angular distribution of quantum scattering is much richer: the features appear even at the angles, where the classical diffraction 
does not exist. N = M = K = 30, lattice period is d = A/2, probe angle is do = 0.17T. 



bution changes depending on K. For example, the distri- 
bution is narrow, when almost all sites are illuminated, 
because the total atom number is fixed and the fluctu- 
ations are suppressed. In contrast, when a smaller part 
is illuminated, the fluctuations are rather large and the 
distribution function is broad. 

If the separated peaks cannot be resolved due to the 
large cavity linewidth, the comb-like structure smooths, 
but its width and form still give the information about 
the atom number distribution function (Fig. 4(b)). In- 
terestingly, in the superfluid state one can get a convolu- 
tion of the Lorentz and Gauss functions, which is the so- 
called Voigt contour, well known in physics of hot atomic 
ensembles with Doppler broadening. Here, it does not 
originates from the motion of hot atoms, but due to the 
quantum uncertainty in the position of ultracold atoms 
[la |. Changing the optical geometry, one can get access 
to the full distribution functions of various atom number- 
related quantities [l6J, e.g., to that of -/V dd — N even . 

Let us summarize the results of this section, which 
established the relations between the properties of the 
quantum sates of atoms and light observables using the 
example of bosonic atoms in optical lattices. The mea- 
surement of scattered light distinguishes between various 
quantum states of ultracold atoms even if their mean den- 
sities are the same. Light scattering is sensitive to the 
global atom fluctuations (i.e. in the whole lattice or its 
significant part), local fluctuations (related to the vari- 
ables at a single site) and distant correlations between 
two or more lattice sites (e.g. the four-point correla- 
tions). Surprisingly, the interesting observation angle is 
not the Bragg one, but the direction of a diffraction min- 



imum, where the strong classical scattering is suppressed 
and the light carries information about the quantum fluc- 
tuations directly. We demonstrated how the Mott insula- 
tor state, superfluid state, and the approximate coherent 
state all can be distinguished by light scattering. The 
use of standing-waves or several probe beams (alterna- 
tively, the quadrature measurements) provides us the in- 
formation about the spatial structure of the atomic noise. 
Moreover, the full atom number distribution functions 
of various atomic variables can be obtained by measur- 
ing spectral (or other phase-dependent) properties of the 
scattered light. 



B. Other models and systems: finite temperatures, 
fermions, spins, molecules 

In the previous section we presented the detailed re- 
sults for the case, where the atomic tunneling between 
neighboring sites was assumed slow, and atomic motion 
was frozen during the measurement. The influence of 
tunneling was explicitly taken into account in the gen- 
eral model in Sec. II. From the general Hamiltonian 
Eq. and Heisenberg equation, Eq. ©, it is evident 
that the dynamics of light and atomic operators are es- 
sentially coupled to each other and the long-range in- 
teraction length can greatly exceed the distance between 
neighboring sites [1JJ . Atomic tunneling affects the phase 
and amplitude of scattered light, and the light scattering 
modifies the tunneling process (e.g. it affects the phase 
and amplitude of the atomic operators). The modifica- 
tion of the light scattering for the time intervals, where 
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FIG. 4: Transmission spectra of a cavity. The spectra directly map out the full atom number distribution functions of an 
ultracold gas. (a) Single Lorentzian for MI reflects the non-fluctuating atom number. Many Lorentzians for SF reflect the 
atom number fluctuations, which are imprinted on the positions of narrow resonances in the spectrum. Here, the cavity is good 
and all satellites are resolved (the cavity decay rate re is smaller than the satellite separation Un, re = O.lUn). N — M — 30, 
K — 15. (b) The same as in (a) but the cavity is worse (re = Uu), which gives smooth broadened contour for SF. Although 
the satellites are not resolved, the spectra for SF and MI states are very different, (c) Spectra for SF with N = M — 70 
and different number of sites illuminated K — 10, 35, 68. The transmission spectra have different forms, since different atom 
distribution functions correspond to different K. K — 0.05[/n. 



tunnelin g is important, was considered, e.g., in 



In Ref. [32[, the short- and long-scale time regimes we 



identified for the case of a BEC in a double-well poten- 
tial. In Ref. [23| , the influence of the photon- induced 
atomic recoil on the site-to-site hopping was analyzed. 
A question of the light dynamics was addressed in Ref. 
34]. 

Here, we consider bosons at zero temperature. A theo- 
retical model of light scattering from bosonic atoms tak- 
ing into account the finite temperature effects was devel- 
oped in Ref. |35| . 

Ultracold gases are not limited to the bosonic atoms 
and models for the interaction of quantum light with ul- 
tracold matter can be generalized for various involved 
many-body states of fcrmions, spin particles, ultra- 
cold molecules, etc. For example, the thermometry of 
fermions in optical lattices based on light scattering was 
suggested in Refs. 36, 37J, where the idea of the light de- 
tection outside the diffraction maxima was shown useful 
as well. 

The consideration of fermions and spins links this di- 
rection of research to the many-body models of con- 



densed matter physics even stronger. Currently, in the 
ultracold gas community, there is a significant interest 
in such many-body phenomena as quantum magnetism, 
fermionic superfluidity, etc. [38| . Quantum light scat- 
tering is a powerful tool to get insight in those effects, 
understand their properties and origins, which still have 
not been fully clarified, and measure their quantum char- 
acteristic. 

Considering the ultracold gas of spin fermions, one 
should underline an important analogy to the system of 
bosons in optical lattices presented above. In the sys- 
tem considered here, the phase of scattered light essen- 
tially depends on the position of an atom and the de- 
tection angle. Such a dependence made possible to ob- 
tain the constructive interference (at Bragg angles) and 
completely destructive interference (at diffraction min- 
ima). More generally, the variation of the observation 
angle was used as a method to tune the relative phase 
of light waves scattered from different atoms. For the 
spin particles, the atomic levels corresponding to differ- 
ent spins can be usually split by some frequency inter- 
val and one can use the spin-selective light scattering 
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for probing those states. Different polarizations and fre- 
quencies of light couple differently to various spin levels. 
The phase difference of light scattered from different spin 
components can be tuned by changing the frequency of 
probing light. Depending on the detuning between the 
probe frequency and two spin levels, on get construc- 
tive or, importantly, destructive interference between the 
light waves scattered from two spin components. Chang- 
ing the probe frequency one can indeed go continuously 
from the totaly constructive to totally destructive inter- 
ference, which is a direct analogy of the diffraction max- 
ima and diffraction minima considered above. One of 
the important consequences of such a method, is that 
at some frequency, the mean light amplitude can vanish 
(destructive interference), and the photon number will 
directly reflect the fluctuations of the population differ- 
ence between two spin levels, thus, revealing the pairing 
between two spin states. In this example, the measure- 
ment of the fluctuations of the atom number difference 
between two levels N^ — JVj. directly corresponds to the 
measurement of the atom number difference between odd 
and even sites 7V 0( jd — N even in a diffraction minimum 
presented earlier in this paper. As the strong classical 
scattering is suppressed in both cases, the light inten- 
sity directly measures the fluctuations of those differences 
given by (iV t - N±) 2 or (iV dd - Ac Vcn ) 2 ■ Creation of spin 
pairs (Cooper pairs) corresponds to the suppressed spin 
difference fluctuations, while for the unpaired states the 
fluctuations are large. As in the case of model presented 
above, this is just one example of the important quan- 
tity, which can be measured. Changing the frequencies 
or angle of light, one can select different variables (as the 
operators Di m ) to be measured in a QND way. 

The pairing of spins into Cooper pairs is a very impor- 
tant phenomenon responsible for the appearance of the 
states such as Bardeen-Cooper-Schricffcr (BCS) or Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) phases. Several the- 
oretical models have been suggested to characterize those 
and other many-body phases of spin particles (including 
quantum magnetism) by light scattering. The propos- 
als involving the QND measurements of spin statistics 
tond the mean-field approximation include Refs. [39l — 
Analogously to the bosons in optical lattices, the 
idea to use the standing wave light for measuring the 
spatial structure of the atomic noise was confirmed in 
Refs. [40, [42[ for spins. As shown in Ref. [40| . various 
quantum antiferromagnetic states could be characterized 
by the polarization-sensitive methods of light scattering. 
While most experimental works on quantum spin fluids 
are currently focused on the mean-field characteristics 
(e.g. using the doubling of the spatial period of spins 
when going from paramagnetic/ferromagnetic to antifer- 
romagnetic phase) [4y, |47| , the experiments (4|| |49| suc- 
cessfully reports the measurements of spin fluctuations 
observed by scattered light. 

Similarly to the case of bosons considered above, the 
cavity transmission spectroscopy can be used to obtain 
the distribution functions of the spin variables in different 



sates of quantum magnets [50H53|. 

Recently, the method of light scattering was applied 
to analyze different systems: ultracold dipolar molecules 
54 5G] trapped in low dimensions (two one-dimensional 
tubes). One finds that in such a system the formation of 
the few-body molecular complexes is possible. In particu- 
lar, the work goes beyond the dimers (similar to those ap- 
pearing in the BCS-BEC theories with two paired spins), 
but also predicts the existence of trimers and tetramers 
(the bound sates of three and four polar molecules, re- 
spectively). The scheme of the light detection at a diffrac- 
tion minimum was generalized and proved effective for 
those few-body complexes. As in the example of bosons 
in an optical lattice considered here, the light scattering 
is sensitive to the position of the particles. It was there- 
fore shown that it is sensitive to the relative distances 
between the molecules forming the complex. The disso- 
ciation (or association) of the complex is accompanied 
by the increase (or decrease, respectively) of the parti- 
cle number fluctuations, which can be readily observed 
in the intensity of scattered light. This opens a way for 
the optical nondestructive in situ detection of the quan- 
tum many-body states of ultracold molecules and their 
dynamics in real time. 



IV. SINGLE-RUN MEASUREMENTS AND 

MEASUREMENT-BASED PREPARATION OF 

THE QUANTUM STATES 

In all problems considered above, the physical mecha- 
nism of light scattering establishes a relation between the 
light and matter observables. As we have shown, the light 
provides the information about various statistical quanti- 
ties of the quantum states of atoms: different correlation 
functions (given by the expectation values of some op- 
erators) and distribution functions of different variables 
can be measured. As it is usual in quantum mechanics, 
the determination of such statistical quantities requires 
multiple measurements. Therefore, the repeated prepa- 
ration of the initial state is necessary, because any quan- 
tum measurement (even a QND one) generally affects the 
quantum state of the system. 

In this section, we will address the problem from an- 
other point of view going beyond the standard goal of 
measuring the expectation values and distribution func- 
tions. Here, we focus on a system dynamics during a 
single run of the optical measurement (i.e. the continu- 
ous detection of scattered photons), without taking the 
average (expectation values) over many realizations as we 
were doing so far [57H60| . Indeed, the result of a single- 
run measurement is important, as it is the first result one 
obtains in an experiment before the averaging procedure. 

During the interaction, the light and matter get en- 
tangled. According to quantum mechanics, due to the 
entanglement, the measurement of one of the quantum 
subsystems (light) will also affect the other quantum sub- 
system (atoms), as an example of the quantum measure- 
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merit back-action. While measurement back-action often 
is an unwanted perturbation, we will use it to change the 
atomic state in a desired way. We will consider the light 
measurement as an active method to prepare particular 
many-body quantum states of the ultracold atoms. 

The Hamiltonian of the system is given by Eq. ([9]). 
Even though some dynamics due the atom tunneling has 
been neglected here, we will show that there is still non- 
trivial dynamics (quantum jumps and non-unitary evolu- 
tion) exclusively associated with the quantum measure- 
ment process. The initial motional state of the ultracold 
atoms trapped in the periodic lattice potential at the 
time moment t = can be represented as 



l* a (0)> = 



E c >i> 



,qu) 



(15) 



which is a quantum superposition of the Fock states cor- 
responding to all possible classical configurations q — 
{<7i, ..,qm} of N atoms at M sites, where qj is the atom 
number at the site j. For each classical configuration q, 
the total atom number is conserved: JZ ■ Qj — A/". This 
superposition displays the uncertainty principle, stating 
that at ultralow temperatures even a single atom can be 
delocalizcd in space, i.e., there is a probability to find an 
atom at any lattice site. We will show, how this atomic 
uncertainty is influenced by the light detection. 

For example, for a limiting case of the MI state, 
where the atom numbers at each lattice site are pre- 
cisely known, only one Fock state will exist in Eq. flT5"l) : 



I^mi) = 1 1, 1 , . -, 1) for the MI with one atom at each 
site. On the other hand, the SF state is given by the 
superposition of all possible classical configurations with 
multinomial coefficients: 



l*SF> 



(VM) 
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E 



AT! 



gi!^!- 
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j\qi,Q2,-qm)- (16) 



Thus, the atom number at a single site as well as the 
atom number at K < M sites are uncertain in the SF 
state. 

We use the open system approach [61| to describe the 
continuous counting of the photons leaking out the cavity. 
When the photon is detected at the moment ti , the quan- 
tum jump occurs, and the state instantaneously changes 
to a new one obtained by applying the cavity photon an- 
nihilation operator \^> c (ti)) — > ai\^c(U)} and renormal- 
ization (the subscript c underlines that we deal with the 
state conditioned on the photocount event). Between the 
photocounts, the system evolves with a non-Hermitian 
Hamiltonian H — ihna\ai. Such an evolution gives a 
quantum trajectory for \^ c {t)) conditioned on the detec- 
tion of photons at times t\, t 2 , .... 

Importantly, due to several assumptions and simplifi- 
cation, one can get a transparent analytical solution of 
the problem. First, after the tunneling is neglected, the 
Hamiltonian does not mix the Fock states in the expres- 
sion (|15[) . So, the problem is reduced to finding solutions 
for each atomic configuration q = {qi, ..,qM}. The full 
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solution will then be given by the superposition of those 
solutions. 

Second, it is known [62, |63[ that, if a coherent 
probe illuminates a prescribed atomic configuration in 
a cavity, the light state is proportional to a coherent 
state \a q (t)) with a q (t) given by Maxwell's equations. 
Each atomic Fock state in Eq. (fT5)) will then be corre- 
lated with a coherent light state with parameters given 
only by the corresponding configuration q: \ty c (t)) = 
J2 q c q exp[^ q (t)}\ qi ,...,q M )\a q (t))/F(t), where F(t) gives 
the normalization. So, the problem to find \^> c (t)) re- 
duces to finding a q (t), $ g (i) for all classical configura- 
tions forming the initial state |^(0)). 

The third mathematical simplification arises due to the 
fact that for coherent states of light, the quantum jumps 
do not lead to a discontinuity of the light amplitude evo- 
lution a q (t) [58J] . Moreover, after t\ 3> 1/k a steady 
state is achieved for all a q (t). Interestingly, in contrast 
to many problems in quantum optics, where the steady 
state is a result of averaging over many quantum trajec- 
tories, here, the steady state in a q (t) appears even at a 
single quantum trajectory. This is a particular property 
of the coherent state of light with damping. 

Although a solution is available for any t (cf. Refs. 
ISMS])' we present it for t > 1/k, when the steady state 
is achieved in all a q (t), and assuming the first photon 
was detected at t\ > 1/k. 

Due to the steady state in all a q (t), the solution is 
independent of the detection times and after m counts is 



induced dynamics of the system, not related to any ob- 
vious processes as tunneling, etc. 

In the following, we will show consequences of Eq. p7|) 
for two cases (I) and (II) addressed in the previous sec- 
tion, where only one probe (ao or 77) exists. As before, 
for the transverse probing (a 7^ 0), we also neglect the 
mode shift, assuming UnD^ <C K or A p . Thus, in both 
examples, a q (|18p depends on the configuration q only 
via a single statistical quantity now called z: z = D q xl 
for cavity probing (77 7^ 0), and z = D\q for transverse 
probing. 

From Eq. (I17[) we can determine the probability distri- 
bution of finding a given z after time t as 



p(z,m,t) 



1 2m -2\a z \ 2 Kt 



Po(z)/F 2 



(20) 
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where the initial distribution po(z) = 
that all configurations q' have the 
12 z \ a z\ 2m Gxp (—2\a z \ 2 Kt)po(z) provides normalization. 
The solution gets especially simple for the case of a 
macroscopic BEC with a large atom number. This situ- 
ation was analyzed in Ref. [59(. 



A. Measurement-induced number squeezing 

As was mentioned, at the Bragg angle, Diq — Nk is 
the operator of atom number at K sites. So, z varies from 
to N reflecting possibilities to find any atom number 
at K sites. The light amplitudes (TT8|) a z = Cz are pro- 
portional to z, C = iUiQao/(iAp — k). The probability 
(j2"0"j) reads 



r] - iU w a D q 10 
iiUnDfi -A p ) + k 



(18) 



p(z,m,t) 



z 2m e- z 



'po(z)/F 2 



(21) 



$q(t) = -\a q \ 2 nt + {rja* q - iU w a Q D q w a q - c.c.)t/2,(19) with a characteristic time r = 2\C\ 2 nt. 



where D q m — ^7=1 u *i{ v j) u m{^j)<lj is a realization of the 
operator Di m at the classical configuration {qi, ..qm}; oo, 
77, and a q all oscillating in steady state at to p were re- 
placed by their constant amplitudes. 

In contrast to a single atomic Fock state, the full so- 
lution, in general, is not factorizable into a product of 
the atomic and light states, and the light and matter 
are entangled. Moreover, in contrast to a single Fock 
state, the quantum jump (applying a\) changes the state, 
and the evolution of the full \^ c (m, t)) is not continuous. 
Note that even for t > 1/k, when all a q reached their 
steady states and are constants, this solution is still time- 
dependent. Thus, the time t — 1/k is not characteristic 
for the steady state of the full solution. 

As we see, each light amplitude a q (t), Eq. (|18[) . is 
given by a Lorentzian corresponding to classical optics. 
Eq. p7[) shows that the probability to find an atom con- 



figuration q, p q (m,t) 



1 2m 



exp(-2K|^)|c^/i^, 



changes in time due to the photodetection. This demon- 
strates the back-action of the light measurement on the 
atomic state and the essentially quantum measurement- 



When time progresses, both m and r increase with 
a probabilistic relation between them. The Quantum 
Monte Carlo method 61] establishes such a relation, thus 
giving a trajectory. Note, that thanks to the simple an- 
alytical solution (fTT|). it gets extremely simple. In each 
step, it consists in the calculation of the photon num- 
ber in the state given by Eq. (|17[) and comparing it with 
a random number generated in advance, thus, deciding 
whether the detection or no-count process has happened. 

If the initial atom number z at K sites is uncertain, 
Po(z) is broad [for the superfluid (SF) it is nearly Gaus- 
sian [16|, |57H59| |]. and Eq. (f2"Tj) shows that p(z,m,t) is 
strongly modified during the measurement. The function 
z 2m exp (— z 2 t) has its maximum at z\ = yrn/r and full 
width at half maximum (FWHM) Sz sa v /21n2/r (for 
Sz <C Z\). Thus, multiplying po(z) by this function will 
shrink the distribution p(z, m, t) to a narrow peak at z\ 
with the width decreasing in time (Fig. 5). 

This describes the projection of the atomic quantum 
state to a final state with squeezed atom number at K 
sites (a Fock states \z%,N — z\) with z\ atoms at K sites 
and N — z\ atoms at M — K sites) . Surprisingly, when 
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Sz < I, the final collapse is even faster than y/r, due to 
the discreteness of p(z,m, t) [601 ■ Measuring the photon 
number m and time t , one can determine z\ of a quantum 
trajectory. 

In contrast to recent results in spin s queezing , which 
can be also obtained for thermal atoms [65l . l66| , in our 
work, quantum nature of ultracold atoms is crucial, as we 
deal with the atom number fluctuations appearing due to 
the derealization of ultracold atoms in space. 

After the distribution shrinks to a single Zi, the light 
collapses to a single coherent state \a Zl ), and the atoms 
and light get disentangled with a factorized state 



|* c ) = \zi,N-zi)\a zl ). 



(22) 



So, light statistics evolves from super-Poissonian to Pois- 
sonian. The conditioned (i.e., at a single trajectory) cav- 



ity photon number (a\ai) c (t) = \C\ 2 Yl z =a z2 P( z > m > *) * s 
given by the second moment of p(z,m,t). Its dynamics 
[very similar to (z) c in Fig. 5(c)] has jumps, even though 
all a z {t) are continuous. In the no-count process, (a\ai) c 
decreases, while at one-count it jumps upwards, which is 
a counter-intuitive signature of super-Poissonian statis- 
tics and conditional probabilities. Finally, it reduces to 



reflecting a direct correspondence be- 



(a\ ai ) c = \Cfzl 

tween the final atom number and cavity photon number, 

which is useful for experiments. 

Importantly, even the final many-body Fock state 
still contains atom-atom entanglement, as many com- 
ponents \qi,..,qM) can have the same z\. For exam- 
ple, the SF state can be represented as \SF)n t M = 
Y, z VbI\ sf )z,k\SF) n ^ Z: m-k {B z are binomial co- 
efficients). After the measurement, it ends up in 
\SF) Zi ^k\SF)n- Z1} m-k, i-e., the product of two uncor- 
rected superfluids. 

Our measurement scheme determines (by squeezing) 
the atom number at a particular lattice region and 
projects the initial atomic state to some subspace. How- 
ever, the atom number at different regions keeps quan- 
tum uncertainty. So, the quantum structure of the final 
state can be revealed in a subsequent optical or matter- 
wave experiment. Thanks to the lattice geometry, one 
can change the illuminated region, and further study 
measurement-induced collapse of the state in the remain- 
ing subspace. 

Even in matter- wave experiments [5j, the product of 
SFs will look different from the initial SF: the atoms from 
different regions will not interfere in average. Note, that 
we did not specify how K sites were selected. One can 
illuminate a continuous region. However, one can illumi- 
nate each second site by choosing the probe wavelength 
twice as lattice period and get number squeezing at odd 
and even sites. In this way, one gets a measurement- 
prepared product of two SFs "loaded" at sites one by 
one (e.g. atoms at odd sites belong to one SF, while at 
even sites to another). While the initial SF shows the 
long-range coherence {b\bj) with the lattice period, the 
prepared state will demonstrate the doubled period in 
(b\bj) (bj is the atom annihilation operator). 



This example shows, that the phase of the matter field 
can be manipulated by this type of QND measurements, 
even though we measure only the occupation number- 
related operators. This is a direct consequence of the 
many-body nature of the atomic state. As usual in the 
QND measurements, by measuring the number-related 
variables, one typically destroys the conjugate variable 
(in this case, the phase). However, here we have a very 
rich choice of the number-related variables for the mea- 
surement. Measuring some of the variables, other vari- 
ables can be left untouched. Therefore, one can carefully 
destroy only particular phase information in our system, 
while other phase-related variables stay unaffected by the 
measurement. This is demonstrated in the last example, 
where the phase coherence survives after the measure- 
ment, but changes its period. 



B. Schrodinger cat state preparation 

In contrast to classical atoms, quantu m g ases scatter 
some light even into diffraction minima [lCl I 111 . Ilq , Il7| . 
Here z = D\ = Y] „_^(— 1) 3+ Qj is the occupation num- 
ber difference between odd and even sites, varying from 
-N to N with a step 2 (we assumed K = M). Eq. (|2"T|) 
keeps its form with a new meaning of z and po(z) [for 
SF, newpo(^) is nearly a Gaussian centered at z = and 
the width y/N 0, M^]. 

The striking difference to the diffraction maximum is 
that our measurement (|21|) is not sensitive to the sign of 
z, while the amplitudes a z — Cz are. So, the final state 
is a macroscopic superposition of two Fock states with 
Z\ t 2 = ±y/m/T and light amplitudes: a Z2 = — a Zl , 

|* c ) = (\ Zl )\a zl ) + (-in -z x )\- a zl ))/V2. (23) 

Figure 6 shows the collapse to a doublet probability 
p(±zi,m,£) and the photon- number trajectory, where 
upward jumps and no-count decreases can be seen. 

In contrast to a diffraction maximum, even in the final 
state, light and matter are not disentangled. However, to 
keep the purity of the state, one should know the number 
of detected photons, because of the sign flip in Eq. (|23|) . 
This reflects the fragility of macroscopic superposition 
states with respect to the decoherence. 

The atom number squeezed states prepared by observ- 
ing light at a diffraction maximum are much more robust 
than the Schrodinger cat state obtained at a diffraction 
minimum to photon lost and detection efficiency, as the 
former do not have any phase jump associated with the 
photocounts. However, the convenient property of the 
measurement at a minimum is that during the same time 
interval the number of photons scattered at a diffrac- 
tion minimum ({a[ai) — |C| 2 iV) is much smaller than 
the one scattered at a maximum ((a[ai) = |C| 2 iV|-) 
[lOL [ill 1571 l58j . Thus the cat state is not rapidly de- 
stroyed by the decoherence as one may expect from the 
coherent Bragg scattering picture. 
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FIG. 6: Photodetection at diffraction minimum leading to the generation of a Schrodinger cat state. Results for a single 
quantum trajectory (Quantum Monte Carlo simulations), (a) Shrinking distribution of the atom-number difference between 
odd and even sites for various times. The doublet corresponds to Schrodinger cat state, (b) Conditional photon number with 
quantum jumps. Initial state: SF, N — 100 atoms, K — M — 100 sites. 



Various types of photon statics during the 
measurement-based state preparation were consid- 
ered in Ref. M. 



C. Quantum state preparation by the 
cavity-transmission measurements 

Let us now consider the prototype configuration (II) 
from the previous section: transmission probing through 
a cavity mirror. In the previous section, we have shown 
that the transmission spectrum has a comb-like structure 
and maps out the full atom number distribution function. 
Such a spectrum is demonstrated in Fig. 4. As discussed 
before, this spectrum reflects the result of multiple mea- 
surements. The question now is what one can get in a 
single run of the optical measurement. Naively, one can 
expect that in each run one gets one of the peaks in Fig. 
4, which would correspond to the collapse to the single- 
peak distribution function. Only after averaging over 
many runs one would recover the full comb- like spectrum 
in Fig. 4. The situation is however much more interest- 
ing in details. Depending on the realization (i.e. for 
different quantum trajectories), the distribution function 
after a single-run of the transmission spectrum measure- 
ment can consist either of one peak or of two peaks. A 
one-peak distribution function corresponds to the atom 
number squeezing, while a two-peak function signals the 
generation of the Schrodinger cat state. 

Thus, probing through a mirror represents an interest- 



ing example, where either a number squeezed state or a 
Schrodinger cat state is prepared depending on the mea- 
surement outcome. Knowing the properties of the system 
after the measurement, one can understand, which kind 
of the states was prepared. Moreover, the cat state ob- 
tained using this configuration has a more general form 
than in the case of the detection at the diffraction mini- 
mum: the amplitudes of two cat components can be non- 
equal and the phase difference between them can take 
various values (not restricted to 7r as in the diffraction 
minimum). The latter property makes the preparation 
of the cat states more robust to the decoherence and 
photon loss [581 ] . 



The full analytical expressions for the probabilities of 
such type of quantum measurement are complicated and 
were reported in Ref. [58]. If for a particular quan- 
tum trajectory the number of photocounts is larger than 
some critical value, the distribution p(z, m, t) collapses to 
a single-peak function centered at z p . This corresponds 
to the collapse to a Fock state \z p , N — z p ) with precisely 
known atom number z p within K lattice sites inside a 
cavity. If, however, the number of photocounts is small, 
the distribution p(z, m, t) collapses to a doublet centered 
at z p with two satellites at Z\2 = z p ± Az. This cor- 
responds to the generation of the Schrodinger cat state 
with generally imbalanced amplitudes and various phase 
shifts between the two components: 
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l*c> = ^[e^+^zuN - Zl )\a Zl )^M^) 

+e-* m *-^\z 2 ,N - z 2 )\a Z2 )y/JM], (24) 

where the phases are tp = — arctan(J7n Az/k) and <&(i) = 
\a Zl \ 2 UnAzt. 

Physically, tuning the probe frequency to A p , we may 
expect maximum scattering for an atom number z p pro- 
viding such a dispersive frequency shift A p = U\\z p . If 
the actual photocount number is large, indeed, the atom 
number is around z p and it collapses to this value. How- 
ever, if m is small, we gain knowledge that the atom 
number z is inconsistent with this choice of A p , but two 
possibilities z < z p or z > z p are indistinguishable. This 
collapses the state to a superposition of two Fock states 
with zi t 2, symmetrically placed around z p . 

The preparation of the particular state is a probabilis- 
tic process. However, importantly, varying the detun- 
ing A p between the probe frequency and the cavity res- 
onance, one can change the distribution of possible out- 
comes and increase or decrease the probability of a cer- 
tain state to appear. In such a way, one can, e.g., increase 
the probability of the cat state with a given imbalance of 
the amplitudes of two components. One can always chose 
the detuning in such a way that the amplitudes will be 
equal to each other independently of the measurement 
outcome. 

In this section, we demonstrated that the cavity en- 
hanced light scattering off an ultracold gas in an optical 
lattice constitutes a quantum measurement with a con- 
trollable form of the measurement back-action. More- 
over, the class of emerging atomic many-body states can 
be chosen via the optical geometry and light frequen- 
cies, as an eigenstate of the operator Di m . The ability 
to implement those operators and, thus, to design the 
measurement back-action is a truly advantage of the op- 
tical lattice geometry in contrast to the simpler setups 
with a BEC in a single- or double- well potential (67H71J , 
where the flexibility is much more restricted. Moreover, 
the use of a cavity selects the prevailing scattering di- 
rection and assures the controllable state collapse, rather 
than complete destroying of the initial atomic state as it 
is the case for the spontaneous emission into all angles 
[72]. Such a measurement-based state preparation is an 
example of the close connection of the measurement and 
decoherence (i.e. the photon leakage and measurement) 
in quantum physics [73, UM ■ 

Considering the projection as a tool for state prepara- 
tion, one can prepare various partitions of the systems. In 
the examples above, one prepared the bipartite systems: 
atoms in the region of K sites and outside that region, 
atoms in odd and even sites. More complicated parti- 
tions are also possible. For example, the light scattering 
from a limited lattice region and detection at a diffraction 
minimum will prepare a tripartite system: the atoms can 
be in odd or even sites, within or outside of the lattice re- 
gion. The initial state of such a system can be expressed 



via the trinomial coefficients. Therefore, using light scat- 
tering, one can simplify the initial problem of numerous 
atoms and sites to the simpler problems of just several 
subsystems. Importantly, the manipulation of particu- 
lar variables can leave other variables untouched by the 
light scattering, allowing the consecutive operations on 
the same sample in the standard QND sense. 

Note, that the quantum state preparation is proba- 
bilistic. However, it can be better controlled by includ- 
ing a feedback loop, which will enable ones the quasi- 
deterministic state preparation. In this case, the trap- 
ping potential should be continuously modified depend- 
ing on the outcome of the photodetector. For exam- 
ple, the detection of photons at a diffraction maximum 
squeezes the atomic number at K sites around some value 
z\, which was not known a priori. The potential can be 
continuously tilted in a way to provide the increase or 
decrease of this atom number to enable ones to obtained 
the number squeezed state with a mean value z\ given 
a priori. The same method can be applied for the atom 
number squeezing at odd or even sites. 



V. QUANTUM OPTICAL LATTICES 

In the previous sections the atoms were assumed to be 
trapped in a prescribed potential formed by stationary 
laser light, which can be described by a c-number func- 
tion. We have shown, how additional light probes can 
measure or modify the properties of the atomic quantum 
state due to the light-matter entanglement. In those ex- 
amples, the quantum natures of both the atomic motion 
and probe light were important, whereas the potential 
was treated classically. 

In this section, we address the more general regime, 
where the quantum nature and fluctuations of the trap- 
ping light potential itself play a key role [12|, uM- Indeed, 
the potential is created by light, which is a quantum ob- 
ject. The goal is to consider the systems and phenomena, 
where the quantum nature of the potential cannot be ne- 
glected and the trapping potential is a quantum dynami- 
cal variable. It must be determined self-consistently with 
the solution for the quantum states of atoms trapped in 
that potential. 

The results of the previous sections can be conceptu- 
ally considered as a perturbation theory for the ultimate 
regime of fully quantum optical lattices. If the trapping 
beams are rather strong, the potential can be represented 
by its mean value and additional quantum fluctuations. 
Then, the quantum fluctuations will obey the equations 
very similar to those for the quantum probe light mode 
considered before. In this way, the entanglement and 
mutual dynamics of light and matter fluctuations can be 
described, including the measurement back-action. Such 
a perturbative regime with small quantum fluctuations of 
the trapping potential can be realized even with optical 
lattices in free space (without a cavity). In this section, 
however, we would like to go beyond such a perturbation 
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FIG. 7: MI to SF phase transition in a quantum optical 
lattice. Cavity influence on the transition is demonstrated 
by the comparison of the occupation probabilities pmi and 
Psf (to find the states |\Pmi) and |*sf}) for a purely quan- 
tum field, i.e., V c i = 0, and a purely classical field, i.e., 
r\\ = 0, as a function of the dimensionless ID on-site in- 
teraction strength g\o I (dEn) (d is the lattice period, Er is 
the recoil energy). We choose rji such that both potentials 
are of equivalent depth, V = 5.5Er, for zero on-site inter- 
action (jid = 0). The quantum (QM) and classical (class) 
cases are depicted with solid and dashed lines, respectively. 
(Un,K,rji) — (— 1, 1/V2, \/5.5)cdR, where the recoil frequency 
is ujr = Er/%. The detuning between the probe and cav- 
ity frequencies (counted from the mean dispersive frequency 
shift) affects the position of the phase transition. In (a) this 
detuning is positive A p — UuN — k and the transition point 
is shifted towards higher interaction strengths in comparison 
to that in a classical lattice. In (b) the detuning is negative 
A p — UuN — —k and the transition point is shifted towards 
the smaller interaction strengths. 



theory and consider the ultracold atoms trapped inside 
a cavity, where the potential is fully quantum. 

A natural implementation to study ultracold quantum 
gas in quantized light is to load the atoms into a high-Q 
cavity. In this case, the light mode of the cavity (e.g., 
a single standing wave, in the simplest case) will form 
the trapping potential. The quantum properties of the 



light mode become more important, when the light in- 
tensity (the number of photons in a cavity) gets smaller 
and smaller. The relative value of quantum fluctuations 
then increase strongly, because the mean light amplitude 
can decrease down to zero. The striking property of 
the cavity configuration is that, although the light in- 
tensity decreases almost to zero, the depth of the optical 
potential in a cavity can be still kept very large. In- 
deed, the potential depth is given by the product of the 
photon number and the light- matter coupling coefficient. 
For example, the quantum potential formed by a single 
standing wave mode a\ with the wave vector k is given 
by hUna\ai cos 2 kx — Ti(g\/ 'A a )a\ai cos 2 kx, where g\ is 
the light-matter coupling coefficient and A a is the fre- 
quency detuning between the light and atomic resonance 
(cf. Eq. §5§ and the discussion after it). In a cavity, 
the light-atom coupling coefficient can reach huge values. 
Thus, even the light field of a single photon can lead to a 
very deep optical potential inside a cavity. The trapping 
and cooling of atoms by a single photon inside a cavity 
has been already demonstrated experimentally [76| (the 
atom in those works was not however cold enough for its 
motion being quantized). 

As we have shown in this paper before, the quantum 
states of atoms and light strongly depend on each other. 
For example, the presence of atoms shifts the cavity res- 
onance and modifies the light amplitude because of light 
scattering. On the other hand, the atoms are trapped 
in that light field, which requires a self-consistent solu- 
tion for the coupled light-matter quantum dynamics. Im- 
portantly, quantum mechanics allows the superpositions 
of several Fock states of photons forming the quantum 
potential. Therefore, one can consider a superposition 
of potentials of several depths. This rises an intriguing 
question about the possibility to have the superposition 
of several atomic phases, each of which is correlated to 
different Fock state of light (i.e. to the potentials of dif- 
ferent depths). For example, one can imagine the poten- 
tial provided by a very week cavity mode, where the Fock 
states of zero photons (the vacuum field), one photon and 
two photons are significant. Then, in principle, one can 
have a superposition, of the free particles (correlated with 
the zero photon number Fock state) , the superfluid state 
(correlated with a potential provided by a single photon) 
and Mott insulator state (correlated with the two-photon 
Fock state). Of course, the stability of such an intrigu- 
ing superposition and its robustness with respect to the 
decoherence should be carefully analyzed. 

Another remarkable property of the quantum optical 
lattices in a cavity is the existence of the long-range in- 
teraction between the atoms beyond the standard Hub- 
bard models. As demonstrated by the Hamiltonian ([7]) 
and Heisenberg equations l[8]). even if the tunneling is 
non-negligible only between the neighboring sites (as it 
is usual for the Bose-Hubbard model), even very distant 
atoms can still interact with each other via the common 
cavity light mode. As a consequence, the tunneling of 
very distant atoms gets correlated (sometimes, called co- 
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tunneling). Mathematically, this is expressed by the co- 
efficients in the generalized Bose-Hubbard model, which 
are not constants, but the operators depending on the 
whole set of light and atomic operators in the whole ex- 
tended optical lattice. Thus, the cavity mediates the 
long-range interaction, whose properties can be tuned by 
the cavity parameters. Even for the simplest type of 
the neighboring interaction, the parameters of the Bose- 
Hubbard model influencing tunneling, atom-atom inter- 
action and atom localization can be controlled by the pa- 
rameters such as probe-cavity detuning, probe intensity, 
probe-atom detuning, cavity linewidth, etc. [12{ . 

As a prominent example, the quantum phase transi- 
tion between the MI and SF states in a fully quantum 
optical lattice was analyzed numerically in a finite sys- 
tem in Ref. [12| • Fig. 7 compares the phase transitions in 
quantum and classical optical lattices of the same mean 
potential depths. From the classical point of view, the 
lattices of the same depth assure the same physics. How- 
ever, in quantum lattices, the potential depth is not the 
only important parameter. In particular, near the phase 
transition point the photon number fluctuations play a 
key role. More precisely, if the transition occurs for the 
potential depth given by the mean photon number n$, 
already the photon numbers n$ ± 1 are associated with 
different atomic phases. Thus the photon fluctuations 
drive the atomic fluctuations and hence the phase transi- 
tion. Depending on the cavity parameters (e.g. the cav- 
ity frequency), the photon fluctuation can either suppress 
or enhance the atomic fluctuations and atomic hopping, 
therefore driving the system towards or outwards of the 
MI or SF states. In Fig. 7 one sees that the position of 
the phase transition in a cavity can be shifted towards 
either smaller or larger values of the atom-atom interac- 
tion strengths, depending on the detuning between the 
probe light and the cavity resonance. (The probe-cavity 
detuning counted from the mean dispersive shift of the 
cavity resonance is positive in Fig 7(a) and negative in 
Fig 7(b)). This demonstrates the tunability of the phase 
transition properties by the cavity parameters. 

Figure 8 compares dynamics of the atomic states in 
classical and quantum lattices of various photon num- 
bers. Importantly, the mean potential depths are chosen 
to be the same, so classically one should not expect any 
difference between all those curves. In contrast, one ob- 
serves a striking difference in the evolutions of the prob- 
ability to find the system in the Mott insulator state for 
different cavity photon numbers. First, for the classi- 
cal potential there is no evolution and the system which 
started in the MI remains in it. For large photon num- 
bers, the MI state gets destroyed with time, but does 
not significantly deviate from the classical solution (be- 
cause the quantum fluctuations are much smaller than 
the mean photon numbers). However, for small photon 
numbers, the fluctuations are larger than the mean value 
and the system leaves the MI state completely. This is 
another example of how the photon fluctuations induce 
the atom hopping and atom fluctuations and radically 
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FIG. 8: Dynamics of the quantum state in a quantum poten- 
tial. Probability pui to find the MI state |^mi) for various 
photon numbers forming the quantum potential. In all cases 
the potential depth is the same (8Er). Classically, the same 
potential depth would mean the same physics. However, in 
a quantum optical lattice, the evolution for the same depth, 
but different photon numbers is strikingly different. For the 
classical potential (constant line) the system stays in the MI 
state without any change. Decreasing the photon number in 
a cavity (the mean photon numbers are 24.2, 4.8, 0.97 and 
0.19), the role of light fluctuations increases, which drives the 
atomic fluctuations destroying the MI state. The smaller the 
mean photon number, the stronger the systems evolution de- 
viates from the case of a classical optical potential [121 ]. 



change the quantum state of the system. Fluctuations 
towards lower photon number enhance tunneling and de- 
stroy the Mott phase. 

The tunneling dynamics in quantum optical lattices 
was analyzed in |13 . |77H73 | . A fundamental difference of 
the quantum lattices from their classical analogues is that 
the coupling of atoms to the leaking light mode opens a 
new dissipation channel for atomic dynamics. In quan- 
tum simulations with classical optical lattices, the main 
advantage of the ultracold atoms in contrast to the con- 
densed matter systems is their almost complete isolation 
from the environment. Using the quantum optical lat- 
tices one can make a step further. The decoherence can 
be introduced in the atomic system in a very careful and 
controllable way: as we have shown in the previous sec- 
tion, the decoherence (the measurement back-action) can 
be tailored using the optical geometry. The results of 
Refs. [H, UM show that the photon number fluctuations 
induce the relaxation of the atom oscillation between the 
sites. More interestingly, looking at a single quantum tra- 
jectory, the tunneling can be completely stopped, when 
the photon is emitted from the cavity, projecting the sys- 
tem to the dark state of the tunneling operator. The next 
photon escape restores the tunneling of the atoms. 

A possibility to generate non-classical states of light 
using the quantum optical lattices at the single-photon 
level was predicted in [801 ] using the few-photon induced 
nonlinearities. 
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The light field in a cavity can show the bistable be- 
havior. Therefore, the trapping potential can be bistable 
as well and different atomic phases can correspond to 
different potential depths. A very exotic phase diagram 
for the MI-SF phase transition in a cavity was obtained 
in Refs. [811, 182). It was shown that the different MI 
insulator lobes in the phase diagram (corresponding to 
the different filling factors) can overlap with each other, 
which is not the case for ultracold atoms in a free space 
optical lattice. The phase diagram significantly depends 
on the cavity parameters. 



VI. VARIETY OF PHENOMENA IN THE 
"ULTRACOLD GAS + OPTICAL CAVITY" 

SYSTEM 

Recent experimental success in trapping a BEC inside 
a cavity [6|-|8j has stimulated a significant attention to 
this system from both experimental and theoretical sides. 
Such a system has already shown a rich behavior for the 
coupled dynamics of the cavity light and atomic motion. 
As not all of those effects require a fully quantum descrip- 
tion and can be explained using mean-field approaches 
either for light or atoms, the detailed review of those 
phenomena is beyond the scope of the present papers. 
However, we will mention some of them here, because a 
deep analysis of the variables, which can carry the in- 
formation about the essentially quantum properties, can 
lead to the determination of particular features available 
only in truly quantum systems (e.g. the light-matter en- 
tanglement). 

One of the interesting aspects of the " BEC + cavity" 
system is its an alog y to the systems known in the field of 
optomechanics [83]. Standard optomechanics considers 
the coup ling between the light and mechanical oscilla- 
tors [8J]- Typically, the setup includes the cavity with 
one movable mirror, where the light mode is coupled to 
the vibrations of that mirror. This simple model was 
extended to various cavities and complicated mechanical 
structures, where the light can be trapped. 

Recent interest in optomechanics was focused on the 
reaching the quantum regime, where the mechanical os- 
cillator should be cooled down to very low tempera- 
tures, which is very demanding experimentally. As it 
was shown, a BEC in a cavity represents a very similar 
system, where the motion of the mechanical oscillator is 
replaced by the motion of the BEC. Indeed, the motion of 
an ultracold atomic ensemble is one of the examples of the 
macroscopic mechanical oscillators. Moreover, the BEC 
is already in its ground state, so the quantum regime for 
such an ultracold oscillator is already reached. Therefore, 
besides representing a significant interest on its own, the 
BEC in a cavity can serve as a system for modeling pro- 
cesses in conventional optomechanics, including its quan- 
tum regime. 

In such novel optomechanical systems, the bistabil- 
ity and nonlinear dynamics have been already demon- 



strated experimentally for the few-photon light in a cav- 
ity [83, Is3- The theories of such effects, including the 
quantum treatment, were developed in J86l - l89j . Mean- 
field approaches for various extensions of that systems 
(e.g. the combined system of a BEC and movable mir- 
ror in the same setup) were considered theoretically in 
90 97]. In the combined system, the bistable behavior 
of the MI to SF phase transition was addressed [98|, [99J . 
The moti on of fermions inside a cavity was analyzed in 
[lOOl llOll ] . Various effects, analogous to the optomechan- 
ical ones can arise in the systems of spins inside a cavity 

doMol. 

A very interesting phenomenon appears, when the 
atoms are illuminated by the light beam transverse to the 
cavity axis. If the light intensity exceeds some threshold 
value, the atoms self-organize in a periodic checkerboard 
pattern occupying only all odd or even sites of the optical 
lattice (the optical lattice is itself a result of the inter- 
ference between the cavity and transverse light modes). 
Such a self-organization phenomenon was proved to be 
linked to the famous Dicke phase transition and the su- 
persolid quantum phase. Usi ng the BE C, this effect was 
obtained experimentally in [1061 . Il07l |. and its p roper- 
ties were analyzed theoretically in Refs. [I08l4ll3| . The 
many- body asp ects of such a system were underlined in 
Refs. Ill 



)d y asp e 

[Mm- 



In this paper we have considered only the far off- 
resonant interaction of light with particles. The main at- 
tention was hence paid to the coupling between the light 
an atomic motion (i.e. the external degrees of freedom). 
Tuning the light closer to the resonance will involve the 
atomic internal excitations. If the ultracold gas is not de- 
stroyed by the spontaneous emission, this will enrich the 
system as at least three quantum oscillators will play a 
significant role. The consideration of atomic excitations 
brings this field closer to the research on e xcitons an d 
polaritons in optical cavities (see, e.g., Refs. [I17l4l2lj ). 



VII. CONCLUSIONS 

Quantum optics of ultracold quantum gases is a quite 
recently developed research direction, which treats light 
and ultracold particle motion at full quantum level. 
Early theoretical proposals, first considered as a kind of 
Gedankenexperiments only, triggered an active and fast 
experimental development, which theories now hardly 
can follow. Particular attention in this field now is paid 
on addressing systems and phenomena, where the use of 
a fully quantum approach is really unavoidable to get 
the physical insight, while the standard mean-field de- 
scriptions of light or matter fail. One should address the 
specific quantities and system properties that carry in- 
formation about the essentially quantum features of the 
physical phenomena. 

Development of this direction can link various fields of 
atomic, molecular, optical and condensed matter physics, 
and even generate the models and approaches for other 
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areas of theoretical and experimental physics. For exam- 
ple, the methods for probing intriguing quantum states 
already discussed can be applied for other systems and 
effects as well : various topolog ical s tates of matter, Dirac 
fermions 112211 . Kondo physics [123J ] , detection o f sta tes in 
fermionic |l24J and superconducting systems [l25J |. etc. 
The methods can be a lso applied for development of laser 
cooling schemes [126l | and precision measurement tech- 
niques |l27l4l3Cll | . The system of quantum gas trapped in 
a cavity provides a significant extension for setups consid- 
ered as candidates for quantum simulators (i.e. atoms in 
classical optical lattices). First, the use of quantum op- 
tical lattices will enable to simulate the broader range of 
trial Hamiltonians. Second, the dissipation can be care- 



fully introduced, tailored and controlled in such a system 
using optical methods. The multipartite entangled states 
of light and many-body atomic system, which naturally 
appear at this ultimate quantum level of the light-matter 
interaction, can serve as a resource for quantum informa- 
tion processing. 
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